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The theory of quasiparticle interference (QPI) for non-centrosymmetric (NCS) superconductors 
with Rashba spin-orbit coupling is developed using T-matrix theory in Born approximation. We 
show that qualitatively new effects in the QPI pattern originate from the Rashba spin-orbit coupling: 
The resulting spin coherence factors lead to a distinct difference of charge- and spin- QPI and to an 
induced spin anisotropy in the latter even for isotropic magnetic impurity scattering. In particular 
a cross - QPI appears describing spin oscillation pattern due to nonmagnetic impurity scattering 
which is directly related to the Rashba vector. We apply our theory to a 2D model for the NCS 
heavy fermion unconventional superconductor CePtaSi and discuss the new QPI features for a gap 
model with accidental node lines due to its composite singlet-triplet nature. 

PACS numbers: 74.20.Rp, 74.55.+V, 74.70.Tx 



The determination of gap symmetry in unconventional 
superconductors (SC) is a persistent problem. Most use- 
ful methods for its investigation are ARPES experiments 
PQ, specific heat and thermal transport measurements 
under rotating field [2] and the use of STM-based quasi- 
particle interference technique (QPI) which utilizes the 
ripples in electronic density generated by random sur- 
face impurities [5H5]- In the cuprate and iron pnictides 
the former is readily applicable but sofar not for heavy 
fermion unconventional superconductors where SC gaps 
are only in the meV range. In this case the last two 
methods are more powerful. QPI technique has recently 
been proposed to discriminate between different d-wave 
pairing states in 115 systems |B] and successfully demon- 
strated for CeCoIn 5 [7] . Before it was also used to inves- 
tigate the hidden order state of URu2Si2 [MTU]. 

Here we propose the application of QPI to inversion 
symmetry-breaking, non-centrosymmetric (NCS) super- 
conductors, notably the tetragonal heavy fermion 131 
and 113 compounds [UJ like CePt 3 Si [12] and CeRhSi 3 
|13j . We develop the QPI theory for the case of mixed 
singlet-triplet gap under the presence of Rashba-type 
spin-orbit coupling. We show that a wealth of new QPI 
features is to be expected: i) Distinct differences in the 
charge- and spin- QPI due to the effect of Rashba coher- 
ence factors, ii) Likewise Rashba - induced anisotropics 
in the spin- QPI even for isotropic magnetic impurity 
scattering, iii) A new kind of cross- QPI where scatter- 
ing by nonmagnetic (charge) impurities leads to spin den- 
sity pattern directly related to the non-zero Rashba vec- 
tor (and vice versa). Finally the expected gap nodes in 
NCS superconductors are generally not on symmetry po- 
sitions but determined by the ratio of singlet and triplet 
amplitudes. QPI can give important information on the 
position of these accidental nodes. We use a weak cou- 
pling BCS theory for the NCS superconductor with a 2D 
model Fermi surface appropriate for 131 compounds. We 
employ Nambu Green's function technique and T-matrix 



theory in Born approximation to calculate the quasipar- 
ticle interference spectra and discuss their new features 
as compared to centrosymmetric unconventional super- 
conductors. 

The BCS model for non-centrosymmetric superconduc- 
tors is given by [TH - fTrj] 
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Here £k = 2ti (cos k x + cos k y ) + 4^2 cos k x cos k y is the ki- 
netic energy with respect to chemical potential /i where t\ 
are next and ti next nearest neighbor hopping. Further- 
more gk = — g-k defines the antisymmetric Rashba term 
due to broken inversion symmetry |17j . Therefore the su- 
perconducting 2x2 gap matrix Ak = [^k^o + dk • cr]ia y 
has an even singlet (?/>k) as well as odd triplet (dk) com- 
ponents. The latter must be aligned with the Rashba 
vector dk = 0ogk to avoid destruction by pairbreaking 
[18] . Here <x = (<j x , a y , a z ) denotes the Pauli matrices. 
Diagonalization of the model leads to an effective two 
band superconductor on the Rashba split bands (£ = ±1) 
given by ekj = £k ~ M + ClSkl that have split Fermi sur- 
faces (FS) and opposite helical spin polarizations as well 
as different superconducting gaps Akj = V'k + £|dk|- For 
concreteness we employ a 2D model for the tetragonal 
Ce- based 131 compounds [TO] [50]. The possible influ- 
ence of magnetic order [5T] is not discussed here. 

In this class gk = g (sin k y , — sin k x , 0) is in the tetrag- 
onal plane. The resulting Rashba-split bands (inset) and 
density of states (DOS) (in the SC state) are shown in 
Fig. [JJl and the electron-type Fermi surface around the 
M(-7r, 7r) points in Fig. [TJj. The split constant-energy sur- 
faces {lu > 0) are presented in Fig. ^p,c. It is known 
from thermal conductivity [22j that the superconducting 
gap has line nodes. To achieve nodes on the M-point 
Fermi surface we use an extended s-wave |18j form for 
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-0k and dk as before. This leads to the two distinct gaps 
Ak£ = ipo + ?pi{cosk x + cos k y ) + £<f> g (sia k% + sinfc^)5 
on the Rashba-split Fermi surfaces = M with quasi- 
particle energies = [ej^ + A?J 2 . The gap zeroes are 
shown as dashed lines in Fig. [ljs-d. Nodes (node lines 
in 3D) appear on the FS for the i^k- quasiparticle sheet 
but not for i?k+ • Their position is accidental, i.e. deter- 
mined by the fine tuning of singlet and triplet amplitudes 
ipi and </>o- The evolution of constant energy surfaces in 
the SC state is shown in Fig.[T]3-h. A few wave vectors q, 
connecting high curvature points close to nodal positions 
that will appear prominently in QPI spectra for small lo 
are shown in Fig.[l^. For larger u (Fig.[TJi) connected FS 
sheets reappear. For the calculation of the QPI pattern 
we use the normal and anomalous matrix Greens func- 
tions of the unperturbed system G = G + a + G_(gk ■ cr) 
and F = [F + ao + ^-(gk ■ c^ify, respectively, where 
(gk = gk/|gk|). 

G±(k ' "" )= 2 £{{}[(»„)» 
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FIG. 1. a) Electronic density of states (DOS) of Rashba 
bands (£ = ±1) in the SC state, (inset: 2D band structure 
along TXMr with t 2 = 0.35t u g = 0.2t 1: and fi = -2t x ). b) 
Normal state electron Fermi surface around M point. Dashed 
lines indicate nodes of gap functions Ak£ with parameters: 
V>o = 2ti, tpi — ti, and <j>o = 0.65ti (also in subsequent fig- 
ures), c) and d) show spectral functions in the normal state 
for different energies oj. e) Spectral function for the super- 
conducting state around M point, are prominent scatter- 
ing vectors defining the QPI patterns (c.f. Figs[3]3, [4^). f-h) 
Spectral function for the superconducting state for different 

UJ. 



-1.5 T8 -4.0 L8 -4.0 L8 -4.3 <L» 




-0.018 0^50 -0.001 0.052 -0.001 0.052 0.001 




FIG. 2. The individual charge- QPI (Ag ) contributions for 
intra- or inter- band scattering with £, £' = ±1, from normal 
(non-magnetic) impurities at the energy ui — 0.2tj: a-d) in the 
normal state e-h) in the superconducting state. (Summation 
over each of the rows is shown in Figplb and FiglSle). 



The QPI density oscillations are obtained from the full 
Green's function that is determined by the effect of scat- 
tering from random charge and spin impurities at the 
surface. The total scattering Hamiltonian in compact 
form reads 

H imp = ^(q)5' Q *k+q/ 5 «*k- (3) 

kqa 

Here we use the Nambu 4-component spinor representa- 
tion ^ = (cj^cJ^c-ktC-H)- In the Nambu space the 
4x4 matrices p a (a — (0,i) = (0,x,y,z)) are given by 
{Pa} = (po,p) = {T3O- ,T a x ,T 3 a y ,T a z ). Here the r Q - 
and a a - Pauli matrices (with t = er = 1) act on Nambu 
and spin indices, respectively. Furthermore we define 
{Sa} = (F S). The first index a = corresponds to non- 
magnetic impurity scattering Vb(q) and entries i = x, y, z 
to isotropic magnetic exchange scattering Vi(q) — V ex (q) 
from impurity spins S. Their components Si are treated 
as frozen, i.e. polarized in a given direction by a small 
magnetic field. An important consequence of the Rashba 
term is that spin and charge channels for impurity scat- 
tering are not decoupled, this also holds true for spin and 
charge response functions |19j . 

The change in STM tunneling conductance in charge 
or spin channel a due to impurity scattering in chan- 
nel (3 is proportional to the resulting change of the elec- 
tron Green's function [3J according to dSI a (r, V)/dV ~ 
-±Im[Tr a p a 5Gi3(r,r,uj = V)] n = 5N a p(r,w). Here 
we focus on the spatial oscillations or momentum de- 
pendence by weak scattering and ignore the possibility 
of bound state formation [23J in the strong scattering 
limit. Therefore we treat the former in Born approxima- 
tion which leads to the Fourier transform of differential 
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conductances given by 
1 



A a/3 (q, iui n ) 



Vs(q)Im A a p(q,iu) n ) 



1 



(4) 
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where N = L 2 is the number of grid points and indices 
(11) refer to the Nambu space. This expression is evalu- 
ated by using Eq. ^ and performing the trace, leads to 
the final result for QPI functions per spin (k' = k q) 



A^„) = ±J2i 1 + ' Sk<)] K$(iun), 

k££' 

AS(^») = E t 1 - ■ 8* - 2gi c gL0]K^(i^ n ), 
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Cgic K ££ q '( iw n), 



(5) 



where the integration kernel for intra-band (£ — (,') an d 
inter-band (£ ^ £') processes is given by 



(iuj n - e^)(iu n - ek-qgQ - A k ^A k _ q g/ 
[(iw n ) 2 -^][(^) 2 -^- qf /] 



These equations give the QPI patterns for the NCS su- 
perconductors. Before presenting numerical results we 
point out the salient new features as compared to cen- 
trosymmetric unconventional superconductors [24] . Here 
A o describes the charge- QPI pattern due to nonmag- 
netic impurities, An the (diagonal) spin- QPI polariza- 
tion pattern due to magnetic impurities and A^o the spin 
polarization generated by nonmagnetic impurities which 
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FIG. 3. a-c) Total charge- QPI (Aq ) in the normal state 
for different cj and scattering from non-magnetic impurities, 
d-f ) The same quantity for the superconducting state. Subset 
(qi,q2,qe) of characteristic QPI wave vectors correspond to 
those in the spectral functions given in Fig[l]e (a; = 0.2ti). 



we term as cross- QPI. It occurs only for finite Rashba 
coupling g k . The coupling of charge and spin degrees 
in cross- QPI appears because the Rashba band states 
are described by hclicity and not spin quantum numbers. 
Therefore in general a particle-hole (charge) excitation 
will also change the spin state. The diagonal QPI func- 
tions in Eq. ^ arc finite even for gk = 0. However 
they are strongly modified when inversion symmetry is 
broken due to the Rashba helicity coherence factors in 
square brackets of Eq. ([5| . The latter have different 
sign for pure charge and spin QPI functions and also 
change sign between intra- and interband contributions. 
The coherence factors can amplify or annihilate the indi- 
vidual contributions in the sum depending on the relative 
orientation of Rashba vectors and band indices. There- 
fore one expects that charge- (A o(q)) and spin- (A^(q)) 
QPI pattern are profoundly different as opposed to non- 
Rashba case where they should be identical. The spin- 
QPI pattern shows a further striking Rashba effect: Al- 
though the exchange scattering itself is isotropic, A^(q) 
is to be expected anisotropic in the present case when gk 
is confined to the tetragonal plane as seen from Eq. jSJ. 
Furthermore the Rashba term leads to nondiagonal ele- 
ments in the spin- QPI pattern, in the present case to 
A^yiiiOn) = A^ x (ioj n ) which is given by 

Kyi^n) =^E^(Slq+SS-,)KgW- (6) 

The QPI functions Aq and A^ are even and A^, is odd in 
q. Therefore 6Nio(r, uj) is finite only when Vb(q) contains 
odd (p-wave) contributions and the net (area integrated) 
cross-QPI conductance always vanishes. 

We conclude that QPI for non-centrosymmetric super- 
conductors derived here exhibits a wealth of new effects 
due to the inversion symmetry breaking Rashba term 
which we present now for the 2D model of 131 com- 
pounds. 

All QPI spectra have four contributions from two in- 
traband and two (equivalent) interband scattering pro- 
cesses. They are shown individually in Fig[2^-d for 
charge-QPI in the normal state and in Fig(2];-h for the SC 
state. In the former the scattering across the Fermi sur- 
face of Fig. [TJd maps out the '2kp<' contours, partly folded 
back into the first BZ. They are slightly different due to 
the different diameter of Rashba split sheets. The diago- 
nal center cross appears only for intraband contributions 
due to scattering process parallel to each Rashba sheet. 
In the SC state for small uj the sheets break up into small 
pieces in BZ regions connecting the node points where the 
gap is small. Therefore the i?k- Rashba band dominates 
because it is the only one with nodes (FigjlJ),e) and the 
SC QPI spectrum reflects the (£,£') = (—1,-1) intra- 
band transitions in Fig(2^. A set of the typical wave vec- 
tors qi of FigjlJ selected by the coherence factors show 
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FIG. 5. a) Total cross- 
QPI of odd spin den- 
sity A^ in the normal 
state from non-magnetic 
scattering, and b) same 
quantity in supercon- 
ducting state, c) Total 
non-diagonal even spin- 
QPI A^ y for the nor- 
mal state from exchange 
scattering by magnetic 
impurities, and d) same 
quantity in supercon- 
ducting state. Here B — > 
is || y and uj = 0.2ti. 



FIG. 4. a-c) Total spin- QPI (A£) for the normal state 
for exchange scattering from magnetic impurities at energy 
uj = 0.2ti. d-f) The same quantity for the superconducting 
state. Each panel corresponds to different orientation of the 
impurity spin S parallel to an applied field B — > 0. Charac- 
teristic QPI wave vectors (Fig[l]e) are shown for B || z. 



up as prominent spots in FigjSjj. The other three con- 
tributions have lower amplitude due to the small £k+ 
sheets at this energy. The sum of all four contributions 
is the total charge-QPI spectrum which is presented in 
Figj3^i-f for the normal and SC state for three energies 
(bias voltages V): uijt\ = 0.1, 0.2, and 0.3. In the for- 
mer (a-c) the '2fcj?' contours which increase with d are 
now blurred due contributions with slightly different di- 
mensions. The diagonal cross due to inraband processes 
appears at all energies. In the SC state for small uj the 
set of prominent spots at q^ due to Fig(2]3 survive and 
are shown explicitly in Fig(3^. The observation of these 
features in QPI and their variation with uj would allow 
to confirm directly the gap model with node structure in 
Fig(T|3,e. For larger uj instead of spots at arcs appear 
due to scattering on and between the again connected 
constant-w surfaces of both Rashba bands in FigjTJi. 

In the centrosymmetric superconductor (gk = 0) 
within Born approximation the charge- and spin- QPI 
functions in Eq. ([5| are identical. The presence of a 
Rashba term introduces two new aspects due to the 
appearance of coherence factors (square brackets in 
Eq. (|5|): i) the charge-and spin-QPI becomes different 
due to different intra-band contributions in the two cases 
ii) While the charge- QPI retains the fourfold symmetry 
as obvious from the first of Eq. ^ and Fig J3j the spin- 
QPI coherence factors in the second Eq. ((5F explicitly 
breaks rotational symmetry with respect to the impu- 
rity moment direction. This is shown in Fig|4] for normal 
and SC phase where we assume that a small magnetic 
field B polarizes the impurity spins S along one of the 
three symmetry directions. The spin- QPI pattern for 



B along x,y directions are still equivalent but rotated 
by 90°, however, both are quite different from the case 
B || z, in particular in the SC state. This is because 
in the present model the Rashba vector fulfils gk • z = 
which singles out the z direction. This anisotropy of spin- 
QPI is therefore characteristic for non-centrosymmetric 
superconductors. The spin-QPI pattern is dominated by 
the full set of scattering vectors q^ connecting the high- 
curvature points of the constant-w surface (Figj4]:) . 

Sofar we have only considered diagonal charge- or spin- 
QPI, however it is a very peculiar feature of NCS Rashba 
systems that charge-spin cross- QPI ( and the reverse) de- 
scribed by Af as well as non-diagonal spin-QPI, e.g. A^ y 
exist. The former describes spin oscillations introduced 
by non-magnetic impurities (and vice versa) which are 
odd under inversion because A^, q = — A^ . It is shown in 
Fig{5^,,b for (xO) component with spin polarization along 
x, again the (yO) case is just rotated in the BZ by 90° 
while it vanishes for (zO). Finally we show in FigJHJ:,d 
the nondiagonal (xy) spin pattern polarized perpendicu- 
lar to the impurity spins and described by Eq. (JgJ) . Only 
the (xy) component is non-zero because gk lies in the 
tetragonal plane. It also exhibits reflection symmetry 
with respect to diagonals (interchange of x,y directions) . 

We have shown that QPI in noncentrosymmetric sys- 
tems exhibits a wealth of new features due to the presence 
of Rashba spin orbit coupling. They are important sig- 
natures of the helical spin texture and a powerful tool 
to investigate the nodal structure of the superconducting 
gap, in particular because the nodal positions are not 
fixed by symmetry for mixed singlet-triplet gaps. 

We thank T. Takimoto for helpful comments and I. 
Eremin for useful discussion. 
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